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ABSTRACT 

A combination of three-dimensional (3D) magnetohydrodynamics (MHD) and syn- 
thetic numerical simulations are presented to follow the evolution of a randomly mag- 
netised plasma that models the intra-cluster medium (ICM), under the isolated ef- 
fects of powerful, light, hypersonic and bipolar Fanaroff-Riley class II (FR II) jets. 
We prescribe the cluster magnetic field (CMF) as a Gaussian random field with a 
Kolmogorov-like energy spectrum. Both the power of the jets and the viewing angle 
that is used for the synthetic Rotation Measure (RM) observations are investigated. 
We find the model radio sources introduce and amplify fiuctuations on the RM sta- 
tistical properties which we analyse as a function of time as well as the viewing angle. 
The average RM and the RM standard deviation are increased by the action of the 
jets. Energetics, RM statistics and magnetic power spectral analysis consistently show 
that the effects also correlate with the jets' power, and that the lightest, fastest jets 
produce the strongest changes in their environment. We see jets distort and amplify 
the CMFs especially near the edges of the lobes and the jets' heads. This process 
leads to a fiattening of the RM structure functions at scales comparable to the source 
size. The edge features we find are similar to ones observed in Hydra A. The results 
show that jet-produced RM enhancements are more apparent in quasars than in radio 
galaxies. Globally, jets tend to enhance the RM standard deviation which may lead 
to overestimations of the CMFs' strength by about 70%. This study means to serve 
as a pathfinder for the SKA, EVLA and LOFAR to follow the evolution of cosmic 
magnetic fields. 
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1 INTRODUCTION 

Combinations of X-ray and radio observations have re- 
vealed the large-scale interactions between radio jets 
from active galactic nuclei (AGN) and the intra-cluster 



llCarilli. Perlev. fc Harris! Il994l : ISmith et all |2002| ) 



dra A IjNulsen et al.l 2002 , and refere nces therein) 



Carilli. Perlev. fc Harris 1991: 



medium (ICM; see e.g. _ 

Blanton. Sarazin. fc McNamaral l2003l : 
2001 l2005l : iNulsen et al.l |2005| ). X-ray surface 



McNamara et al.l 



bright- 
ness depressions, or X-ray cavities, are seen in the 
ICM, the position of which correlates with that of 
the synchrotron emitting radio galaxy lobes at GHz- 
frequ en cies, e.g. NGC 1275 in Per se us ijBoehringer et al.l 
1 19931 : ISchmidt. Fabian, fc Sanders! l200aV Cygnus A 



Hy- 

and 

several other clusters (jBirzan et al. I l2008l ). AGN inject large 
amounts of energy into the ICM in the form of relativistic, 
magnetised and coUimated plasma jet s. These seem to be 
activ e for periods of some tens o f Myr (^Alexan der fc Leahy! 

jets 
. itjy 

heated llBinnev fc TaboJ !l995! : !Kaiser fc Alexander. .19991: 



^ e t ens o i Myr f 

! 19871 : !Kaiser fc Alexandei! [l999l : !Krau sc 2005). The jets 
drive st rong shocks in the I C M which is co nsequently 



Reyn olds, Heinz fc Begelman 2001; Alexandeil 2002! : Krausd 
!2003h . Jet plasma and magnetic fields inflate a cavity, 
the "cocoon" , w hich expands and displaces the ambient 
medium plasma (!Scheuerl 1974) ; Kaiser fc Alexandeil !l997! : 
!Heinz. Reynolds, fc Begelman 19981 ). The material inside 
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cocoons is responsible the observed radio synchrotron 
emission. 

The interaction between AGN jets and the ICM has 
been m odelled as a hydrody namical process jClarke et alj 



19971: iKaiser fc Alexander 'l999'; 'Vicent et alJ I2OOII: 



Churazov et alJ [2OOI.: Alexander .2002.: Brtiggen fc Kaiseil 



2OO2I : iBasson fc Alexandej l2003l : iKrausd 



200 



It is well 

established however that magnetic fields thread the ICM. 
Radio polarimetry is the major source of information about 
cluster magnetic fields (CMFs). T he AGN radio emission 
is typically ~ 10-50 % polarized (|Bridle fc PerlevI 1 19841 ) 
and the i ntervening magnetised ICM gives rise to Faraday 
rotation (|Kim. Kronbe rg. fc Tribble' 1991, and references 
therein; for a review see Fcretti & Giovannini 2008) , which is 
characterized by the Rotation Measure (RM). Maps of this 
quantity have been produced using radio sources located 
within the ICM or be hind it, as well as usin g cluster radio 
halos (for a review see ICarilh fc T avlor"2002Y The distribu- 
tion in RM images, along with the assumption that CMFs 
vary over a single scale, has lead to estimates of the strength 
of CMFs within 5-30 ^G in the co ol-core clusters Cygnus A 



Hydra A, A1795 and 3C 295 fseelDreher Carilli. fc Perlev 
19871 : iTavlor fc PerlevI 1 19931 : ICe fc Owenll l993l : I Allen et al 



200 ll . respectively) , within 2-8 /x G in the non-coo l core 



clusters Coina, A 119, 3C129 fsee iFeretti et al.lll995l . Il999l : 
iTavlor et aLll2q0ll. res pectivelvl and also in A2634 and A4Q0 
ijEilek fc Owenll2002l ). The RM also suggests that CMFs 
located in the cores of clusters are stronger than those 
in their outer regions (iKim. Kronberg. fc Tribbld Il99ll : 



IClarke. Kronberg. Bohringer 200 ll ). A correlation between 
the RM values and the mass deposition rates of cluster 
cooling flows has been detected as well (|Soker fc Sarazii] 
Il990l ). 

The structure of CMFs seems to be mainly related to 
turbulence. This has been inferred using the spatial varia- 
tion of observ ed F araday rotation images and their statisti- 
cal properties (fmbb le 1991: Enfilin fc Vogt 200 |), in partic- 
ular, the shape of the RM s tructure funct i ons |Laing et al.l 
|200^ : iGuidetti et al] I2OI0I '). iMurgia et al. used the 

RM images of A119 and modelled its magnetic fields with a 
three-dimensional multi-scale tangled geometry. They used 
spatial variations of the field characterized by a magnetic 
power spectrum following a power law in Fourier space. Mur- 
gia et al. found that both the mean and the standard de- 
viation of the RM, which are used to estimate the strength 
of CMFs, depend on the magnetic spectral index, n. They 
conclude that n ~ 2 produces the bes t fit to the observa - 
tions in A119. With a similar approach ICovoni et al.l (|2006l ) 
reported spectral power law slope fits of n ~ 2 at the core 
of A2255, and of n ~ 4 at the outer par ts. Also, indexes clos e 
to 11/3 have been estimated in A2382 (Guidctti ct al. 2008"), 
A400, A2634 and Hydra A (|Vogt fc Enl31in.2003 , 2005 ). 

Models of the evolution of CMFs suggest that the fields' 
structure is shaped by the dynamics of the ICM via magnetic 
fiux freezing an d ideal MHD which is a good approximation 
in clusters fsee |jone3l2008l , and references therein). In this 
context, CMFs should follow the ICM advection closely, 
particularly the one driven by AGN jets. Turbulence in 
the ICM is expected to be produced during cluster for- 
matio n and mergers (iRuzmaikin. Sokolov. fc Shukurovl 



19891: iRoettigerrStone. fc Burns! 

Dolag. Bartelmann. fc Lesch 20021 : Schindleil 



1999 
2002 



[Dubois fc Tevssiej I2OO8I: IRvu et all l2008h. as weU as 
by A GN outfiows (jEnfilin fc Vogtl I2OO6I : iRuszkowski et all 
|2007| . and references therein). Turbulent ICM fiows are 
thus believed to amplify CMFs from weak seed fields to 
the observed microGauss r egime by the formation o f local 
small-scale dynamos (see ISchekochihin et al.l |2004| . for a 
numerical review). 

Characteristic depolarization gradients are often ob- 
served between the ma in hotspots o f seve r al po werful 
radio galaxies (see e.g . iLaingj Il988l ). iLaingj |l988) and 
iGarrington et ahl l|l988l ) have explained such features using 
path length differences between each of the radio lobes and 
the edge of the observed X-ray emitting ICM that contains 
them. The RM map o f the Fanaroff-Riley Class I (FR I, 
iFanaroff fc Rilevlll974h radio source Hydra A shows oppo- 
site signs in the appro a ching and the rece ding radio lobes 
(|Tavlor fc PerlevI 1 19931) . iLaing et al.l l|2008D have described 
these RM fluctuations with synthetic three-dimensional 
Faraday screens characterized by (i) a random distribution 
of magnetic fields, (ii) a spherically symmetric gas distribu- 
tion, (iii) axisymmetric prolate-spheroid cavities resembling 
the geometry of the Faraday screen projected by Hydra A. 
These 3D differential Faraday screen models are consistent 
with the radio source expansion framework of ,Scheuer. (.1973 ) 
and also with observations of X-ray cavities. The RM gradi- 
ents in question, therefore, imply that AGN jets make their 
way by pushing the ICM gas and its magnetic fields a t once . 
This process was recently studied bv^Guidetti et al] (|201ll ) 
with high resolution multi-frequency images of Faraday rota- 
tion and depolarization of the radio galaxies 0206-1-35, M 84, 
3C 270 and 3C 353. The latter is an extended Fanaroff Ri- 
ley Class II (FR II) radio source such as the ones that we 
model in this paper. Guidetti et al. reported, for the first 
time, well defined RM bands, with very little small-scale 
structure, perpendicular to the major axis of the radio lobes 
of these galaxies. Also, they use a simple time-independent 
model of both the radio source compression and the ambi- 
ent medium magnetic fields to generate the RM bands in 
synthetic observations. 

Observations of CMFs are limited by the resolution and 
sensitivity of teles copes. Many questions about these fields 
are still open (see lBeck. Brunetti. fc Feretti. L.II2006I . for a 
review). At least for powerful FR II sources there seems to be 
little mixing of ICM plasma with the synchrotron-emitting 
one; the observed RM distri bution is caused by plasma ex- 
ternal to the r adio source ( Dreher. Carilli. fc Perlevlll987l : 
ICarilli fc Taylor 2002). The interaction of expanding radio 
sources with the magnetised ICM should happen locally 
to the sources. An important question here is how signif- 
icant this process is when interpreting RM observations 
as indicators of the properties of both the CMFs and the 
ICM gas themselve s. The SKA (jHuarte-Espinosa et aLll2009l : 
iKrause et al1l2009l ). EVLA, e-MERLIN, LOFAR and XEUS 
(the new generation of telescopes) are planned to have ultra- 
high sensitivity and sub-arc-second resolution which will al- 
low a better understanding about the cosmological evolu- 
tion and structure of CMFs. e.g. we will unders tand impor- 
tant details about the ICM turbulence (see e.g. iLaing et al.l 
2008^; 'Guidctti et al."2010") and its heat conducti on (see e.g. 
Schckochihin ct al. 2005; Bogdan ovTc et al.ll2009l ). 

This is the second of two papers in which we use nu- 
merical simulations to study the evolution of magnetic fields 
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in Fanaroff-Riley class II jets and their immediate cluster 
surroundings. In this work we introduce a combination of 
3DMHD and synthetic RM numerical simulations. We fol- 
low the large-scale interaction of powerful radio jets from an 
FR II radio source with the magnetic fields in the core of a 
non-cool core cluster. The structure, the evolution and the 
polarised synchrotron emission of the radio sources them- 
selves are discussed in our earlier paper (Huarte-Espinosa, 
Krause & Alexander, 2011; Paper 1). Here we concentrate 
on the structure of the ICM and the feedback from radio 
sources. We only follow the evolution of the sources while 
jets are active. The behaviour of both energetics and mag- 
netic power spectra (in Fourier space) of CMFs is explored 
as a function of position in the two-dimensional parameter 
space defined by the velocity and density of the jets. Syn- 
thetic RM maps are produced at different viewing angles to 
investigate the observational signatures of expanding radio 
sources. 

This paper is organised as follows: in Section[2]we briefiy 
describe the formalism of ideal MHD, the numerical methods 
we use and write about our implementation of the ICM, the 
CMFs and AGN jets. The results of our simulations are pre- 
sented in Section[3]aIong with their interpretation and anal- 
ysis from the point of view of energetics, RM statistics and 
magnetic power spectra in Fourier space. Both the implica- 
tions and applications of our work are discussed in Section U 
There, we also examine the well-studied case of Hydra A and 
finish by commenting about our model assumptions. Finally, 
in Section [S] we summarise the main results of our work and 
present our conclusions. 



2 SIMULATIONS 

2.1 Governing equations 

To describe the dynamics of the plasma in the ICM and AGN 
radio jets, we use the system of nonlinear time-dependent 
hyperbolic equations of ideal compressible MHD. In three di- 
mensions and non-dimensional conservative form, these are 
given by: 

g+V(pV) = ft (1) 

^^ + V- (pVV + p + bV2-BB) = Pg + Pj(2) 

l^+V- [(iJ + p + BV2)V-B(V.B)] = is, (3) 

^-Vx(VxB) = 0, (4) 

where p, p, V and B are the plasma density, thermal pres- 
sure, fiow velocity and magnetic fields, respectively. In Q, 
E = p/ipi - 1) + pV^ /2 + and is the total energy den- 

sity, while 7 is the ratio of specific heats. In the right hand 
side of (O, ((2]) and ([3]) we have source terms to implement 
jets via the injection of mass, pj, momentum, Pj, and kinetic 
energy, Ej (see Section \2A\ . as well as a Newtonian gravi- 
tational acceleration, g, in order to keep the initial plasma 
in a magneto-hydrostatic equilibrium. 



Table 1. Simulations and parameters. 



Simulation 










name 


[Mach] 




[xlO^s W] 


[Myr] 


lighter-slow 


40 


0.004 


4.6 


14.1 


light-slow 


40 


0.020 


17.2 


8.3 


lighter-fast 


80 


0.004 


28.1 


7.1 


light-fast 


80 


0.020 


128.8 


4.4 


lighter- faster 


130 


0.004 


112.8 


4.7 


no-jots 








all of 










the above 



Time average average jets velocity in the nozzle. It is equal 
to the external Mach number. 

^ Time average average jets vs ambient density contrast in 
the nozzle. 

^ Jets power from equation (11) in Paper 1. 
Simulations' ending times. 

2.2 Code and implementation 

We solve the equations of ideal MHD in three dimension s 
using the numerical code Flash 3.1 (|Frvxell et al.l 120001 ). 
Flash's new multidimensional unsplit constrained transport 
(CT) solver is employed to m aintain th e divergence of mag- 
netic fields down to < 10"^^ jLee fc De anc 2008). A diffu- 
sive HLLC solver (|Batten et al.ll 19971 ) prevents spurious low 
pressure and density values from appearing in the grid. We 
use a Courant-Friedrichs-Lewy parameter of 0.25 and peri- 
odic boundary conditions in all the domain's faces. These 
boundary conditions prevent numerical noise from polluting 
the magnetic spectrum in the grid (Section l2.3.ip . Our com- 
putational domain is a cube with edges |x| ^1/2 computa- 
tional units (which represent a volume of 200 kpc"^), Carte- 
sian coordinates and a uniform grid with 200"' cells. 

We present six simulations to study the effects of ac- 
tive AGN jets on the magnetic fields in the core of a non- 
cool core cluster. Five simulations, which we refer to as jets- 
simulations, are designed to experiment with the power of 
the jets. We do this by varying their velocities and densi- 
ties. These jet-simulations are the same that we present in 
Paper 1. The other simulation, the no-jets model, is car- 
ried out in order to follow the evolution of the ICM plasma 
without any perturbations throughout the same simulation 
timesteps as the jet-simulations. All of these models are la- 
belled descriptively and summarised in Table 1. See Paper 1 
for a discussion about the low densities and hypersonic ve- 
locities of the jets. 

2.3 Initial conditions 

We implement the cluster gas using an ide al gas equation 
of state, a King density profile (|Kinel[l972l ) and a state of 
magneto-hydrostatic equilibrium with a background Newto- 
nian central gravity field (see Paper 1 for details). 

2.3.1 Cluster magnetic fields 

The magnetic eld within the cluster is set up as an isotropic 
random eld with a power law energy spectrum Following 
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iTribbld lIlQQll ) and iMurgia et al.l l|2004 ), we generate a cu- 
bic grid in Fourier space. For each cell, we define three 
components of a vector potential which takes the form 
A(k) = A(k)e''''''', where k is the frequency vector (fc^ = 
k^ + k^ + kl) and i is the unitary complex number, while A 
and 9 are the vector's amplitudes and phases, respectively. 
We draw 9(k) from a uniform random distribution within 
and 27r, and A(k) is also randomly distributed but has a 
Rayleigh probability distribution 



P{A,e)dAd6 



A 



2-K\Ak\ 



-exp 



A^ 



dAd6, 



where we choose the power law Ansatz 



(5) 



(6) 



for a given slope C,. 

We transform to real space by taking the inverse fast 
Fourier Transform jPress et al.lll993 ) of A(k). The resulting 
magnetic vector potential A(x) is multiplied by the plasma 
density radial profile (see Paper 1, Section 2.3.1). This multi- 
plication implements the magnetic flux freezing by generat- 
ing magnetic fields, the strength of which follows the plasma 
density and pressure profiles. 

We then take the curl of this vector to get the magnetic 
field which is normalised so that the ICM thermal pressure is 
approximately ten times larger than its magnetic pressure. 



B2/2 



10, 



(7) 



which is a reasonable value in this context (|Carilli &: TavloJ 
I2OO2D . 

This procedure yields solenoidal magnetic fields which 
are tangled at scales of order our computational resolution, 
and characterized by spatial variations that follow a mag- 
netic power spectrum with a power law of the form 



-C-t-2 



= k" 



(8) 



We choose a Kolmogorov-like three- dimensional turbulent 
slope n = -11/3, based on the w ork of lVogt fc EnfilinI l|2003l . 
[iooi) and lGuidetti et all (|2008l '). Note that other values for 
n would also have been possible (compare Section 1, above). 
This choice may also be motivated by theoretical stud- 
ies o f the cosmological evolution of cluster magnetic fields 



fe.g.lBriiggen et al 



ISchekochihin et al.l 



200, 



200. 



Rvu et al.ll20 08: but compare also 



Xu et al.ll2010l '). The topology of 



these fields is (of course) constrained by the approximation 
of ideal MHD and thus we do not model physical magnetic 
reconnection, but the simulations do present numerical resis- 
tivity. We note that the Fourier method implicitly imposes 
maximum and minimum scales. 



2.4 Jets 

By implementing source terms to the equations 
([2)1 and ([3} we inject mass, momentum and kinetic energy to 
the central grid cells located within a control cylinder of ra- 
dius, rj, and height, hj, which are resolved by 3 and 8 cells, 
respectively. Inside this "nozzle" we update the plasma den- 
sity and x-velocity using the constant source terms pj and 

Vy 

Bipolar back-to-back jets are continuously injected until 



they reach the computational boundaries. The simulations 
are stopped at this point. The nozzle plasma pressure, pj, 
takes the constant value of the central ambient pressure (i.e. 
Pc/7), while HyVyjV^ and E evolve according to the MHD 
equations Q - @ . This implementation is similar to that of 
lOmma et al.l |2004i ) . but generalised here to the MHD case. 

The power of the jets is given in Table 1 and computed 
using the equation (11) in Paper 1. 

2.4-1 Cocoon contact surface 

We implement a passive tracer, r(x,f), which is injected 
along with the jets' plasma to distinguish it from the one 
of the ambient medium. Defined as a real number within 
(Ix 10^^'',.99), r initially takes these extreme values in the 
ICM and in the jets' nozzle, respectively. The tracer is ad- 
vected with the jets, and a comparison of the distribution of 
both r and p allows us to identify the contact surface of the 
cocoon with an accuracy of about 4 computational cells. 



3 RESULTS AND ANALYSIS 

To analyse our simulations in the context of galaxy cluster 
cores, we scale the plasma variables in the following way. The 
size of the computational domain is given by L = 200 kpc, 
the ICM speed of sound as Cs = 1.66x10^ kms~^ -which 
corresponds to a temperature of 10* K - and the clus- 
ter central density as pc =1 x 10~^^ kgm~^. Each compu- 
tational cell is one kiloparsec long, the central cluster ra- 
dius is = 160 kpc and the nozzle's radius and height 
take the values Vj = 3 kpc and hj = 8 kpc, respectively. For 
the given ICM properties, the initial RMS magnetic field 
is ~ 20 microGauss. 



3.1 Evolution of the plasma 

We let the magnetised ICM plasma described in Section 2.3 
to relax for ~ 118 Myr before we inject the jets. This 
timescale represents one thermal crossing time (i.e. the ra- 
tio of the computational domain's length over the gas sound 
speed) . During this phase the magnetic energy of the cluster 
dominates the kinetic energy which was initially zero. CMFs 
produce local, subsonic and homogeneous random motion on 
the plasma. We see the three-dimensional magnetic power 
spectrum of the plasma preserves the initial Kolmovogov 
turbulent profile. Mild losses of about half an order of mag- 
nitude develop in the spectrum at scales < 10 kpc due to 
numerical difi'usion. We only follow this phase in order to 
reduce the initial differences between the magnetic and the 
thermal pressure profiles of our initial setup (Section 2.3), 
but we do not intend to simulate the evolution of the ICM 
and its magnetic fields at any time before the injection of 
our jets. 

The left column in Figure 1 shows the ICM plasma at 
the end of the relaxation phase. These are logarithmic false 
colour maps of the plasma density (top row), the magnetic 
field strength (middle row) and the magnetic beta parameter 
(bottom row). The injection of the jets starts at this point; 
at tjet =0Myr. Then we see that our simulations produce 
the basic fluid dyn amical structures found in previous jet 
simulations (e.g. see lBasson fc Alexandeij2003l : [Omma et al.l 
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[a} logipicsi ) , *j = 1 = Myi 




(d) lo(}{\B\^)Ajet ^OMyr 





(b) log<j>£CM), M£ich=40, r\ —0,02 & 
tjed =a,3Myr 





(c) logifiicM), Mfi<±i=13£i, ij =D.0O4 Sc 
t,,i =4.7Myr 




(f) roi?(|J3p), Madi-^1,W, jf r-O.OCM & 
t_,,e -4,7Myr 




(i) io5(^)„J,Mach= 130, fj =0.004 & 



Figure 1. Logarithmic false colour maps of the plasma density (top row), the magnetic fields strength (middle row) and the magnetic 
beta parameter (bottom row). The top and middle rows are shown in computational units. The left column shows the ambient plasma 
after one relaxation time (equivalent to 118.5 Myr), but just before the injection of jets. The middle column shows the plasma after 
the effects of the light-slow source at the end of the simulation. The right column shows the ambient medium after the injection of the 
lighter-faster source at the end of the simulation. 
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Figure 2. Central slice thought our cubic computational domain 
showing the geometry of cluster magnetic field lines under the 
effects of the "lighter-fast" source at tjot =5.3Myr. The cocoon 
is marked with yellow dots and jets are in the plane of the image, 
horizontally. We note the grayscale has no relation to the strength 
of the fields at all. 



I2OO4I : iRevnolds. Heinz, fc BegelmanllioOll : lKrausdl2005h . A 

strong bow shock is formed in the ambient medium, behind 
which a contact discontinuity bounds the cocoon inflated by 
the jets' former plasma and magnetic fields. The middle col- 
umn in Figure 1 shows the plasma after the effects of the 
light-slow source at the end of that simulation. In the right 
column, we present the plasma after the injection of the 
lighter-faster source at the end of that simulation (see Ta- 
ble 1). Note the dependence of the cocoon geometry, which 
is marked by yellow surfaces in Figures 1(b) and 1(c), on the 
density of the jets. A full discussion of the source properties 
and evolution is covered in Paper 1. 

As the sources develop, the ambient plasma is com- 
pressed and the shocked region between the cocoons and the 
bow shocks expands. During this process the CMF compo- 
nent that is normal to the direction of the cocoon expansion 
is compressed and stretched (Figure [2]). The strength of the 
CMFs is thus amplified, particularly ahead of the jets' work- 
ing surfaces, and shows local enhancements above the mean 
field intensity -Figures 1(e) and 1(f). Such enhancements 
reach factors up to about 7.4, 7.7, 10.9, 7.5 and 12.2, re- 
spectively, for the jet-simulations in the order in which they 
appear in Table 1. 

We see (Figure 1, bottom row) that jets produce consid- 
erable gradients in Pm inside the sources; the cocoons and 
the ambient medium show values of order 10000 and 10, 
respectively. Yet the global change of fim. in the ambient 
medium is fairly low, Pm ^ 10 throughout the simulations. 



3.2 Energetics 

We have calculated the energetics of the ICM as a func- 
tion of the source size which is measured from hotspot to 
hotspot. Figure [3] has three panels which show the ther- 
mal, the kinetic and the magnetic energies of the ICM, from 
top to bottom. Each panel has six profiles corresponding to 
the simulations in Table 1. We see that in the absence of 
jets (black dotted curve) the ICM thermal energy increases 
slightly in time due to dissipation of magnetic energy which 
falls by ~ 10% due to numerical diffusion. When jets are 
present, on the other hand, the total energy of the ICM 
increases consistently. The energy transfer from jets to the 
ambient medium is predominantly thermal and proportional 
to their velocity. Jets with a density contrast of 4xl0~^ are 
more efficient in this process than those with a density con- 
trast of 2xl0~^ (Figure [HJ. This behaviour is characteristic 
of very light jets which inflate fat cocoons that provide large 
areas over which to affect the surrounding ICM. In contrast, 
the forward momentum of relatively heavier jets leads to 
smaller cocoons. Also, both the pressure and the speed of 
sound inside the cocoons with lighter jets are high er than 
those inside the cocoons with relatively heavier jets ijKrausd 
l2003h . 

The decay of the ICM magnetic energy is slowed down, 
but not suppressed, by the Mach-40 jets (see the green and 
dark-blue curves in Figure |3}. In contrast, the faster jets are 
able to signiflcantly impede the decay by the time that their 
cocoons extend for ~70kpc (from hotspot to hotspot). As 
the lightest, faster sources continue expanding (see the red 
and light-blue curves) they increase the energy of the CMFs 
in proportion to the jet velocity. This increase however tends 
to flatten modestly after the cocoons of these sources are 
130 kpc long. 



3.3 Rotation measure 

We use the information from our simulations to produced 
synthetic RM maps and compare the predictions of our mod- 
els with the observations. Our maps are produced by com- 
puting 

D/kpc 

RM^812 / (£^)(^)'^'-d/m^ (9) 

Do/kpc 

along the line of sight, where Do,D,ne and are the lo- 
cation of the computational domain's boundary, the loca- 
tion of the cocoon's contact surface (the spatial distribution 
of which is given by the tracer of the jets; Section l2.4.ip . 
the ICM electron density and the CMF component that is 
parallel to the line of sight, respectively. We note that the 
magnetic fields inside the sources are not considered for this 
analysis, and hence, there is no internal rotation. 

The viewing or inclination angle, 9v, is measured from 
the jets' axis to the line of sight. Jets point towards the ob- 
server when 9v = 0° , and higher viewing angles are obtained 
by rotating the simulation data cubes perpendicularly to the 
jets, anticlockwise. Thus jets are on the plane of the sky, 
horizontally, at 9^ = 90° . For a consistent comparison at all 
viewing angles, the computational cells which are located 
at radii larger than 100 kpc are discarded from our cubic 
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Figure 3. ICM energetics as a function of source size measured 
from liotspot to hotspot. 



computational domain of 200 kpc'' (Figure |4]). The geome- 
try of the Faraday screens, or the RM integration volumes, 
depends on both the shape of the cocoons' contact surface - 
which is different in each simulation - and the viewing angle. 
Observationally, the measured screens are modified by noise 



and the finite telescope resolution as well (jEnfilin &: Vogtl 
l2003h . 



3.3.1 No-source RM maps 

We are only interested in the statistical changes of the RM 
which are caused by the expansion of FR II sources. In our 
simulations however we also see modest, yet detectable, RM 
variations that are caused by subsonic random flow in the 
ICM, independently of the source expansion (see top of Sec- 
tion [3]T}. In order to interpret the mechanical work specif- 
ically done by the sources on the ICM's gas and magnetic 
fields, we produce "no-source" RM maps. To this end we use 
equation (|9]), but combine information from different simu- 
lations in the following way: 

• He and B\\ are taken from the no-jets simulation (Sec- 
tion g^J, 

• Do and D, which define the geometry of the Faraday 
Screen, are taken from the jets-simulations (see Table 1). 

This idea is illustrated in Figure [l] Hence "no-source" RM 
maps give the RM structure which would have been observed 
against an evolving radio source had it not affected the ICM 
dynamically. 

We have produced no-source RM maps at tjct = 0, 4.4, 
4.7, 7.1, 8.3 and 14.1 Myr (i.e. just before the injection of 
jets and at the five timesteps that correspond to the end of 
the jets-simulation, te; see Table 1). 



3.3.2 Synthetic RM maps 

Figures [51 4101 show our synthetic RM maps and their his- 
tograms. Each Figure corresponds to one of the jets- 
simulations and a viewing angle (either 45 or 90 degrees). 
Panels are arranged by rows. The top row shows (from left 
to right): the no-source RM map at tjot =OMyr (using the 
cocoon coordinates of the corresponding source at tjct ~ to); 
the no-source RM map at tjct ~ to (also using the cocoon co- 
ordinates of the corresponding source at tjet ~ to); the RM 
map at tjot = to- The middle row, (b), has panels in which 
the black, the red and the green profiles correspond to hori- 
zontal lines through the left, the middle and the right maps 
in the top row, respectively. The middle panel in row (b) 
shows 1-D RM profiles at j/ = kpc. The left and right pan- 
els in (b) show 1-D RM profiles at symmetric j/-values with 
respect to y =Okpc, which go near the sources' edge. Fi- 
nally, the bottom row, (c), shows histograms of the three 
maps in the top row - where (again) the black, the red and 
the green profiles in (c) correspond to the left, the middle 
and the right maps in (a), respectively. 

Our synthetic RM maps consistently show that pow- 
erful jets enhance the RM distribution locally. We see this 
process scales with both the jet velocity and density. More 
spherical cocoons, such as those that are typically shown by 
very light jets, generally yield higher RM values; they favour 
the alignment of enhanced magnetic field components with 
the line of sight. 

The considered RM enhancement occurs because the 
expanding hypersonic sources compress the ICM gas and 
magnetic fields in the shocked ambient region which is lo- 
cated just ahead of the cocoons' contact surface. Magnetic 
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Figure 4. Sketch of the geometry employed for the no-source RM maps. The Unear colour scale shows the y component of the CMFs 
normalized to their maximum value in the X-Z plane through the centre of the cubic computational domain. No-source RM images are 
produced using JSJ; gas density and magnetic field factors are taken from the no-jets simulation, whereas the integration limits are taken 
from a jets-simulation. See Section 3.3.1 for details. 



field lines are compressed and stretched perpendicularly to 
the source expansion direction (Figure [2]). We see that this 
RM increase is typically larger at the sources' edge, where 
the alignment between the line of sight and the compressed 
cluster magnetic field lines is more important. Moreover, 
shocks are stronger near the hotspots than near the sides of 
the cocoons. Such structure is also evident in the RM maps 
-e.g. see the "lighter- faster" source map, Figure |8] (a), right 
panel. 

By comparing the top left and top middle panels in 
Figures [S tilOl as well as the black and the red profiles in 
the rows (b) and (c), we see a modest RM variation in the 
"no-source" RM maps. Relative to such small changes, the 
source produced RM enhancements are larger, particularly 
those of the jets with velocities exceeding Mach 80. The "no- 
source" RM maps produced at 6^ ~ 45° show the expected 
global asymmetry; having a shorter total Faraday depth, the 
approaching source lobes consistently show lower RM values 
and variation than the receding lobes. We see powerful radio 
jets enhance these RM gradients and that the histograms in 
panel (c) of Figures [5t 410l show that this effect happens at 
the highest RM values. When jets are on the plane of the sky 
[dv = 90°) the RM histograms tend to be double peaked. Yet 
at 9v = 45°, the histograms are narrower and single peaked; 
the Faraday depth is more uniform at 9v ~ 90°. Further, 
the alignment of the compressed CMF components with the 
line of sight is more important for observations at angles 
within 20-70 degrees due to the prolate spheroid geometry 
of cocoons. Therefore, we see higher RM values at 45 degrees 
than at 90. The choice of 45 degrees is interesting because 



it separa tes quasars fro m radio galaxies in the unification 
model of lBarthej l|l989h . 



3.3.3 RM evolution 

The time evolution of the RM statistics has been followed 
as the sources expand. In Figure[TT]the top row displays the 
evolution of the mean RM. The bottom row shows the evolu- 
tion of the RM standard deviation. The viewing angle is 45 
and 90 degrees for the left and right columns, respectively. 
RM maps are produced using expanding sources at corre- 
sponding timesteps. Profiles are normalised to their corre- 
sponding initial values (i.e. to either the mean RM or the 
RM standard deviation produced against Faraday screens of 
matching sources, but on the gas and magnetic field distri- 
bution at tjct ~0, see Section [3. 3.ip . 

As the sources expand we see that the RM statistics 
scale up with the speed of the jets. Generally the "lighter" 
sources yield more important changes on the statistics than 
the "light" sources, by factors of about 2-7. We see the gradi- 
ents in the RM statistics evolution profiles (Figure [TT|) follow 
similar trends to those of the ICM energetics profiles (Fig- 
ure Also, in agreement with the RM histograms (panel 
c, Figures [5l410p , the effect of the sources on the RM evolu- 
tion appears to be more significant for observations made at 
45 degrees than at 90. This is, again (see previous section), 
caused by the geometry of cocoons. 

Gradients in the RM statistics evolution profiles (Fig- 
ure [11]) also show a relation with those in the pressure evo- 
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(a) No-source RM map at t =OMyr (left), No-source RM map at t =7.1Myr (middle), RM map at t =7.1Myr (right). 
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(b) RM at y = —25 (left), y =0 (middle) and y =25 kpc (right). The black, red and green curves correspond to the left, the middle 
and the right panels in (a), respectively. 




(c) RM histograms of the three maps in (a). Colors are 
as in panel (b). 



Figure 5. RM maps and histograms of the lighter-fast simulation at 9^ = 90°. (a) From left to right; the no-source RM map at 
tjot = OMyr (using the cocoon coordinates of the lighter-fast source at tjot = 7.1 Myr); the no-source RM map at tjet = 7.1 Myr; the RM 
map at tj^t = 7.1 Myr. (b) Each panel shows RM profiles corresponding to horizontal lines through the three maps in (a). 
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(b) RM at y = —25 (left), y =0 (middle) and y =25 kpc (right). The black, red and green curves correspond to the left, the middle 
and the right panels in (a), respectively. 
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Figure 6. Same as Figure [5] but for 9y = 45°. The approaching lobe is on the positive distance range. 
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Figure 7. Same as Figure [5l but for the lighter- faster jets. 
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Figure 8. Same as Figure[5]but for the lighter-faster jets and 9v = 45°. The approaching lobe is on the positive distance range. 
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(b) RM at y = —12 (left), y =0 (middle) and y =12 kpc (right). The black, red and green curves correspond to the left, the middle 
and the right panels in (a), respectively. 
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Figure 9. Same as Figure|5]but for the light-fast jets. 
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(b) RM at y = —12 (left), y =0 (middle) and y =12 kpc (right). The black, red and green curves correspond to the left, the middle 
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Figure 10. Same as Figure[5]but for the light-fast jets and S„ = 45°. The approaching lobe is on the positive distance range. 
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lution profiles that we shown in Figure 1121 The latter show 
that the "slow" sources reach pressure equilibrium with the 
ambient medium earlier than faster sources {vj ^ 80 Mach) . 
In agreement we see that the "slow" sources cause modest 
changes on the RM, relative to faster sources. We see that 
the "lighter-faster" source, which has the most powerful jets 
amongst our models, enhances both the mean RM and the 
RM standard deviation by approximately 70% for a viewing 
angle of 45 degrees (red curve. Figure [TTl left column). This 
effect is very important because, as we will discuss in the 
next Section, the observed values of the RM dispersion are 
used to infer the strength of CMFs. 

3.4 Field structure evolution 

In order to investigate the effects that jets have on the 
structure of CMFs, we calculate three-dimensional magnetic 
power spectra in Fourier space using the output from our 
simulations. Specifically, we use the magnetic fields of the 
five jets-simulations which were produced a.t t — (see 
Table 1). The spectra are shown in Figure 13(a) For this 



computation we exclude the volume that is occupied by the 
sources (we use the spatial distribution of the jets' tracer. 
Section [2XT|) . 

In general we find that the spectra preserve the initial 
condition of a power-law with a Kolmogorov index (Sec- 
tion [2]3TT]). Yet the energy injection from the jets produces 
power enhancements on scales that range from the resolu- 
tion limit, 1 kpc, to the source size, ~40kpc. An interesting 
feature that we see is a flattening on the power spectra at the 
scale defined by the cocoon size - Figure |13(a)| This effect 
depends on the power of the jets. As radio sources expand, 
the energy of the jets is transferred to the ICM gas and 
CMFs in the shocked ambient region. There, magnetic field 
lines are compressed and stretched tangentially to the sur- 
face area of the radio cocoons (Section l3.ip . Hence, we find 
that fatter or more spherical cocoons, such as the ones pro- 
duced by the "lighter" jets, are more efficient in increasing 
and redistributing the energy of the CMFs. 

The scales at which the jets affect the CMFs are more 
evident in Figure 13(b) There we plot the ratio of the power 
3(a) over the corresponding spectra com- 
no-source" simulation, i.e. for each plot 



spectra in Figure 
puted using the 



in Figure 13(b) we use the corresponding profile in Fig- 
ure 13(a) as the numerator. The spectrum in the denom- 
inator is calculated in the following way: (i) we combine 
the magnetic fields of the "no-source" run at tjct = te with 
the cocoon spatial distribution of the corresponding jets- 
simulation at tjct ~ te- (ii) The magnetic fields inside the 
cocoon are removed. This yields a cube filled with magnetic 
fields, which have not been affected by the jets at all, and 
an empty cavity in the middle. The geometry of such cavity 
exactly matches the shape and orientation of the cocoon's 
contact surface, (iii) We compute the 3D magnetic power 
spectrum of these fields. We note the numerical diffusion of 
the simulations is cancelled in the computation of the pro- 



files shown in Figure 13(b) 



Gradients in these spectral ratio profiles persistently 
show maxima close to 40 kpc, as well as a roughly constant 
and progressive increase from scales close to 10 kpc down 
to Ikpc -see Figure 13(b) These features are related to 
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Figure 12. Evolution of the average cocoon thermal pressure as 
a function of the source size measured from hotspot to hotspot. 
The ICM thermal pressure is also shown (bottom dashed black 
line) for reference. 



to the thickness of the ambient medium swept-up region. 
This gas component expands as the sources evolve towards 
a thermal pressure equilibrium state with the ICM gas (Fig- 
ure [T2ll, and reaches scales of about 10 kpc at the end of the 
simulations. 

We see that the expansion and size of our model co- 
coons are clearly linked to the enhancement of the CMF 
power spectra at large-scales. The cocoons which are formed 
by the lightest, fast jets are bigger and have more isotropic 
gas and energy distributions (see Section |3]2]). Figure[2]illus- 
trates the deformation of the CMFs' topology that is caused 
by the "lighter-fast" source at 9^ = 90° and tjot = 5.3 Myr. 
The cocoon is marked by dots and it clearly stretches and 
compresses the neighbouring magnetic field lines, particu- 
lar ly_alMad_o£_th«jets^_^ surface. 

'En filin fc Vogtl { 20031 ) have demonstrated that there is 
a correlation between the statistical properties of the RM 
distribution, in the real space, and the structure of magnetic 
fields in the three-dimensional frequency space. Based on 
these facts we compute the RM structure functions of those 
of our synthetic RM maps that correspond to tjct = te and 
9y = 45° . The structure function is given by 



5(A) 



7V(A) ^ 

(x ,y) .{x' ,y' ) 



[RM{x,y)-RM{x',y')]\ (10) 



where the sum is taken over the A''(A) non-zero RM pix- 
els that are inside bins (an nuli) centred at (x, y) and have 
radii of A = \{x,y) — {x',y')\ (jSimonetti. Cordes. fc Spangled 
Il984h . 



Figure 13(c) shows the RM structure functions in ques- 



the source energy and momentum injection scales and also 



tion. At the end of the simulations, the structure functions 
generally still show the power law that we have used as the 
initial conditions of the simulations fSection l2.3.ip . Yet these 
curves considerably flatten at scales of order tens of kpc. The 
structure functions rise again from a certain scale on. We 
see that this scale is larger for that sources with fatter co- 
coons. The flat regions therefore seem to be directly related 
to the cocoon width. This is consistent with our findings 
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Figure 11. Evolution of the RM statistics as radio sources expand. Panels at the top row display the evolution of the mean RM, whereas 
panels in the bottom show the evolution of the RM standard deviation. Panels in either the left or the right column correspond to synthetic 
observations made at 45 or 90 degrees, respectively. The RM is computed using expanding sources at corresponding timestcps. Profiles 
are normalised to the mean RM, or the RM standard deviation, of corresponding no-source RM maps at tj^t =0 (see Section l3.3.1l l. 



Table 2. Summary of the RM statistics evolution. These are the maxima of the mean RM 
and the RM standard deviation, normalized to their corresponding initial values, reached in 
synthetic observations that capture the effects of powerful radio sources on CMFs, in terms 
of the jet power and 6v ■ 



Simulation 


max((/?Af> / (-RAfini)) 


max((T/(Tini) 


max((_RM>/(M/ini)) 
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(b) Ratio of the power spectra in (a) over the power spectra of 
corresponding no-jcts simulations (see Section [3.4l for details). 
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4 DISCUSSION 

A combination of 3D-MHD numerical simulations and syn- 
thetic RM observations have been produced and presented. 
These simulations follow the evolution and observational sig- 
natures of CMFs under the effects of powerful jets from 
an AGN in the core of a non-cool core cluster. To under- 
stand the basic physics of this interaction, the effects from 
other cluster gala^xies, synchrotron cooling, magnetic recon- 
nections and observational factors (such as beaming, depo- 
larization and light travel effects) have not been considered 
in the model. We have focussed on the role that the jet ve- 
locity, the jet density and the observational viewing angle 
play on the jets-CMFs interaction. Plasma energetics, RM 
statistics and magnetic power spectra consistently indicate 
that CMFs are affected by AGN jets. The effect depends 
on the jets' properties; the lighter and the faster the jets, 
the more important the results they have on the ambient 
medium. 

The shape of the CMF energetics profiles (Figure [3| 
is consistent with the following arguments. Radio sources 
expand due to the energy injection from jets. The cocoons 
displace the ICM gas and magnetic fields from the cavity, 
compressing them into the region located between the con- 
tact surface and the bow shock (Figures [T] and [2]) . These 
CMF lines are also stretched in this process. The magnetic 
component tangential to the surface of the radio cocoons 
is therefore enhanced. As the sources expand, the cocoon 
pressure decays and the compression caused by bow shocks 
weakens (Figure [T^ . This is in good agreement with radi o 
source evolution models (e.g. see lKaiser fc Alexandeilll997l ). 
Once the sources with Mach 40 jets (green and dark-blue 
curves. Figure (T^J have expanded for about 100 kpc, their 
cocoon thermal pressure is comparable with that in the en- 
vironment. Thus only these relatively slow sources are not 
powerful enough to significantly affect the magnetic energy 
of the ICM. 

Our RM integration neglects the contribution from ICM 
gas and magnetic fields beyond the computational domain. 
Jet-produced RM enhancements are local, important at the 
edge of the sources and should not depend on this outer 
RM component. Observations show that both the clus- 
ter ion density and the s trength of the CMFs decreases 



away from the cluster core jKim. Kronberg. fc Tr ibblc 1991 
IClarke. Kronberg. BohringeJ l200ll : ICarilh fc"Tavlor. .2003 . 



5/3 



and references therein). Thus we expect to capture the major 
contribution to the RM by our procedure. Since the region 
we use for the analysis is spherical, any additional contri- 
bution to the RM is expected to be statistically spatially 
uniform. 



10""! I 

1 10 100 

Separation [ kpc ] 

(c) RM structure function produced at 8v = 45° 

Figure 13. Structure of CMFs. (a) and (b) are in the three- 
dimensional frequency space, (c) is in the two-dimensional scale 
space of the RM. 
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4.1 Characteristic RM gradients 

iGuidetti et al.l (|201ll ) have recently reported bands in the 
observed RM images of the FR II radio sources 3C 353 and 
Cygnus A, as well as in the FRI sources 0206+35, M84 
and 3C 270. These bands are perpendicular to the major 
axis of the radio lobes and are characterized by very little 
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Figure 14. Profiles of RM statistical fluctuation projected onto 
the jet axis of our "lighter-faster" source. The viewing angle cor- 
responding to these maps is 45° . The receding lobe is in the neg- 
ative distance range where fluctuations arc the strongest. 20 kpc 
represent 18.7arcsec approximately. 



small-scale RM structure. To explain the nature of such fea- 
tures, Guidetti et al. present synthetic RM observations gen- 
erated with a simple model of gas with constant density and 
magnetic field distributions. This plasma is then affected 
by compression induced by an elliptical cavity that is artifi- 
cially placed inside the gas. The model allows these authors 
to conclude that there are three magnetic field components 
associated with observed RM: an isotropic field which has 
large-scale fiuctuations, plausibly associated with the undis- 
turbed ICM; a well-ordered field draped around the leading 
edges of the radio lobes; a field with small-scale fiuctuations 
which is located in the shells of compressed gas surround- 
ing the inner lobes. We do not use an ordered magnetic 
field component in our model and do not observe RM bands 
in our synthetic images. Thus we confirm the findings of 
Guidetti et al., namely that without any ordered field com- 
ponent in the ICM, the radio source expansion will not lead 
to banded RM structures. 

In addition, our model shows that source-induced RM 
variations can be as large as 70%, are correlated with the 
radio jet properties, depend on time and are stronger at the 
edge of the cocoons. This is a very important point because 
the observed RM signal is used to infer the strength of the 
CMFs (e.g. seelFeretti et al.., 19951 . Il999l : iTavlor et al"]|200ll : 
lEilek fc Qwenll2002l ). 



4-1.1 Asymmetrical RM gradients between radio lobes 

Our synthetic RM maps persistently show pronounced 
gradients at regions corresponding to the receding ra- 
dio lobes (the left ones) relative to the approaching ones 
(Figures [s UlOp . Such RM differences have been observed 
in the FR lis Cygnys A JPreher. Garilh. fc PerlevI ligSTl- 



Carilh. Perlev. fc Harris! Il994l ) and 3G 353 l|Guidetti et all 
as well as in several FR Is (e.g. 3G 31 and Hydra A, 
Laing et all l2008l: 3C449. iGuidetti et all I2OI0I : 0206-1-53, 
3C 270 and M 84, iGuidetti et al.ll201ll '). Although we simu- 
late RM enhancements caused by FR II radio sources, FR Is 
-which are known to form ICM X-ray cavities as well- 
may well produce similar effects. Hydra A is a well studied 
central cluster radio source. The Faraday rotation in this 
case is thought to be dominated by the local ICM contribu- 



tion. This cluster presents an ICM magnetic spectral index 
of n Rill/3 which has been inferred against the emission 
from the embedded luminous FR I radio source Hydra A. 
The inclination angle of this source has been inferred to be 
within 37-60 de grees (T aylor fc Perlevlll993l : lEnfilin fc Vog^ 
I2OO6I : iLaing et al... 2008 1 . 

In terms of RM statistical fluctuations that are asym- 
metrical between the source's lobes, we find that our 
"lighter-faster" simulation ( Table 1) agrees to a c ertain 
degree with the studies of iTavlor fc PerlevI |l993') and 
iLaing et al.l l|2008l ) on Hydra A. This is partly unexpected, 
as we have chosen a particular realisation for the magnetic 
power spectrum of the CMFs. However, we focus here on the 
edge effects, which is where we believe that radio sources 
have the strongest effect: Two panels are shown in Fig- 
ure [14] where we present curves of (RM) (left panel) and 
CRM (right panel) as a function of position along a hori- 
zontal line between the lobes. The mean and standard de- 
viation of the RM are taken on vertical slices of the map 
corresponding to the "lighter-faster" simulation at 6^ = 45° 
and tjot =4.7 Myr. Vertical error bars are computed by tak- 
ing the difference between the red and the black curves in 
Figure[8] Both statistical measures have large gradients near 
the left edge. This effect is produced by the expanding radio 
source. Towards the right edge, the radio source expansion 
has also strongly amplified the RM values. However, the 
initial magnetic field in that region happened to be small, 
which is why this region does not look spectacular. This 
is quite comparable t o the RM di s tribut ion in observed ra- 
dio sources, e.g. see iLaing et al. 1 l|2008l . Figure 22, panels 
a and b) who present th e RM observations of Hydra A 
made by ITavlor fc PerlevI (119931 ) along with synthetic RM 
observation. It is therefore well possible that such large RM- 
gradients towards the edges of radio sources are strongly in- 
fiuenced by the lobe expansion. Note that we consider here 
only one realisation for the random Faraday screen. A de- 
tailed match of the considered RM gradient with any obser- 
vation is therefore not intended. 

Our model jets inject energy at scales that range from 
the resolution limit to the source size, but preferentially 
near the latter. As we have seen, this results in a fiattening 
in both the magnetic po wer spectra and th e RM structure 
functi ons at large-scales. Laing et all (|2008l ). [Guidetti et al.l 
l|2010l) and IGuidetti et al.1 (|201ll ) have reported such flat 
gradients in the observed structure functions of the FR I 
sources 3C 31, 3C270, 3C 449 and Hydra A. Our simula- 
tions show that the physics of radio source expansion may 
explain these structure function gradients naturally. Thus 
conclusions about the nature of the ICM turbulence may 
not be drawn from the gradients in question. 



4.1.2 The Laing-Garrington effect 

We have calculated the mean and the standard deviation 
of the RM for the approaching {{RM)'^'^ and cr'^) and the 
receding {{RAiy and o-^m) source lobes, as a function of 
the viewing angle. In the top left panel of Figure [TJ] we show 
(i?A/)''P in sohd lines and {RMy in dashed lines. In all of 
these plots, positive (negative) angles correspond to inte- 
grations along (against) the line of sight, and again jets are 
on the plane of the sky at ±90°. In the bottom left panel 
of Figure [15] we show o-^m/<^iim- addition, to demon- 
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Figure 15. RM statistics comparison between tiie approacliing and receding radio lobes as a function of tlie angle between 
the jet axis and the line of sight. Top: mean RM of the approaching (solid) and the receding (dashed) lobes. Bottom: 
""iiM (receding lobe) /(Tfljv/(a'PP'"oaching lobe). Positive and negative angles represent RM integrations along and against the line of 
sight, respectively. Panels in the right column show data from no-source RM maps (Section 13.3. 11 1. 



strate the effect of the jets on the Faraday screens, which 
will in turn affect the RM computations, we produce anal- 
ogous {{RMy, {RM)''^ and ct^'m/'^rm) plots but use the 
data from the no-source RM maps (Section l3.3.1|l . We show 
these no-source RM plots in the right column of Figure 1151 
side by side to their (with-) source counterparts. 

Together, all of these profiles show that the general de- 
pendence of the lobe RM statistics on the viewing angle is 
shaped by the intrinsic nature of the Faraday screens, as 
it is generally assumed in the interpretation of the Laing- 
Garrington effect. The result of the radio source expansion 
is to amplify several RM features. 

Only for positive inclinations our realisation results in 
stronger average RMs in the receding lobes (Figure 1151 
top panel). We see strong amplification at inclinati ons of 
20° < < 70°. Q-RM is linked to the depolarisation (|Burnl 
ll966l : lKrause et al.llioOTi '). Already intrinsically, the depolar- 
isation is always higher for the receding lobe (o-rm/ct^'J^^ >1). 
Panels in the bottom row (Figure [TSj) show that this trend is 
mostly amplified by the radio sources, apart from the fatter 



(lighter) ones for which the cr-ratio drops below unity near 
&v = ±90° . Thus the Laing-Garrington effect is only moder- 
ately affected by the radio source expansion, in such a way 
that the associated trends tend to be amplified. 

We have shown that powerful radio sources affect the 
structure of CMFs which is characterized by the magnetic 
power spectrum index n. This process is important because 
the ICM energy transformation, ei ther from or to mag netic 
form, as well as the observations (iMurgia et al.ll2004l ). de- 
pend on n. Ifence we agree with iLaing et al. ( 20081 ) in that 
the structure of CMFs inferred by using the emission of clus- 
ter radio sources should take a piecewise defined form over 
the range of the spatial frequencies. Our simulations show, 
moreover, that such functional forms change with time (see 
Figure 113b . b; note that the profiles in panel b would be 
completely flat otherwise). Faraday rotation measurements 
from clusters with broad differences in age will be key in 
the determination of evolving magnetic structure functions 
of the CMFs. 
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5 SUMMARY AND CONCLUSIONS 

We have presented a combination of 3D-MHD numerical 
simulations and synthetic RM observations to model the 
evolution and observational signatures of cluster magnetic 
fields under the effects of powerful AGN jets at the core 
of a non-cool core galaxy cluster. We prescribe the cluster 
magnetic field as a Gaussian random field with a power law 
energy spectrum tuned to resemble the expectations for tur- 
bulence. We focus on how both the jets' power, in terms of 
their density and velocity, and the observational angle affect 
the observed RM signal. Such signal is used by observational 
studies to infer the strength of cluster magnetic fields. Our 
results are very consistent and show the following. 

Powerful jets increase the ICM energy in proportion 
to their velocity and in inverse proportion to their density. 
Light jets are more efficient in this process because they 
inflate fatter cocoons than relatively heaver jets. Only jets 
with Mach numbers of 80 and 130, which are over-pressured 
with respect to the ambient medium, significantly increase 
the ICM magnetic energy. We see that the time evolution 
of the CMF energetics is correlated with that of the cocoon 
pressure. 

Jets compress and displace the magnetised ICM. This 
results in a flattening of the 3D magnetic power spectra and 
the associated fleld RM structure functions. This happens 
at scales from the resolution limit to the source size, but 
preferentially near the latter. We see the effect is correlated 
with the jet velocity and pronounced for fat cocoons. The 
RM structure functions of 3C 31, 3C270, 3C 449 and Hy- 
dra A show flat gradients at large scales, and the cause of 
such shapes has been prev iously modeled (|Laing et al.ll2008l : 
iGuidetti et al]|2010l . [20T1I ) in terms of ICM turbulence. Our 
model suggests, however, that the physics of radio source 
expansion may as well explain this feature. 

The general distribution of the Faraday depth is de- 
termined by the intrinsic nature of the ICM. Powerful jets 
enhance the RM locally by factors up to ^1.73. We flnd 
this process is proportional to the jet velocity and in in- 
versely proportional to the beam density. Expanding hyper- 
sonic sources compress the ICM gas and magnetic fields in 
the shocked ambient region. The strength of the fields that 
are perpendicular to the local lobe expansion direction in- 
creases, specially near the hotspots. The alignment of com- 
pressed CMF components and the line of sight is favoured 
(i) at the cocoon edge (ii) by sources with a fat cocoon (iii) 
at 9-u within 20-70 degrees, due to the cocoon geometry. 
Jet-produced RM enhancements are thus more important 
under these conditions and should be more evident in radio 
galaxies than in quasars. 

Both (RM) and aRM continuously change and tend to 
increase as the radio sources develop. "Lighter" sources with 
jet velocities ^ 40 Mach yield {RM) and o-rm enhancement 
of factors within 1.2 to 1.7, in proportion to the jet velocity. 
These RM enhancements may therefore lead to overestima- 
tions of the CMF strength by up to about 70%. 
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